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ABSTRACT 

The suppression of synchrotron emission at low frequencies 
due to the influence of the ionized medium is investigated. 
Explicit expressions and detailed numerical values of the 
emission spectra are presented for various electron energies 
and plasma and cyclotron frequencies. Unlike previous studies 
of this suppression effect, which were applicable only to 
ultrarelativistic electrons, the present treatment is valid 
for electrons of arbitrary energies and in particular for 
intermediate energy electrons such as those presumably ac- 
celerated in solar flares. 
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INTRODUCTION 

The suppression of the low frequencies of synchrotron 
emission due to the influence of the ambient ionized medium 
was investigated by Tsytovich (1951) , Ginzburg (1953) and Razin 
(1957, 1960). This suppression effect is essentially the result 
of the phase velocity of light becoming greater than c at frequencies 
above the plasma frequency but sufficiently near to it, where the 
index of refraction of the ambient medium is less than unity. Ad- 
ditional studies of this effect and its application to cosmic 
radio emission were given by Ginzburg and Syrovatskii (1964, 1965), 
Scheuer (1965) and McCray (1966, 1967). 

The influence of the ambient coronal plasma on the spectra 
of solar Type IV radio bursts, which are generally believed to 
be synchrotron emission of energetic electrons accelerated in 
solar flares, was investigated by Ramaty and Lingenfelter (1967). 

They showed that the low frequency cutoffs observed for these 
bursts ( Takakura and Kai, 1961) could result from the suppression 
of the synchrotron emissivity at low frequencies rather than from 
absorption during propagation through the solar corona, and that 
from the study of the observed cutoff it is possible to deduce 
the magnetic field or the ambient density of the emitting region. 

This possibility was recently substantiated by Boischot and Clavelier 
(1967) who pointed out that the very sharp low frequency cutoff 
of a Type IV solar burst, which they observed to originate at a 
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large distance from the solar surface, could indeed be a sup- 
pression effect rather than absorption, since at the frequency 
of the observed cutoff, the ambient coronal plasma could hardly 
affect the propagation of radio waves. From these considerations, 
they then provided an estimate of the coronal magnetic field at 
the site of the emission. 

The straight-forward application of this suppression effect 
to solar radio emission, however, must be treated with some caution 
because all the available treatments of the influence of the ionized 
medium on synchrotron emission are valid only for ultrarelat ivist ic 
electrons whereas the bulk of the Type IV bursts could originate 
from electrons of only mildly relativistic energies (Takakura, 1960; 
Holt and Cline, 1968). 

The radiation in vacuum of electrons of arbitrary energies 
spiraling about the lines of force of a constant magnetic field 
was first derived by Schott (1912) and discussions of his treat- 
ment were given by Schwinger (1949) and Landau and Lif shitz (1962) . 
Takakura (1960) used these results to calculate the synchrotron 
emission from intermediate energy solar electrons. His formalism 
was applied to various Type IV bursts (Takakura and Kai, 1961; 
Takakura , 1967; Holt and Cline , 1968), but, as was pointed out 
recent 1 y ^Takakura , 1967; Holt and Cline , 1968) this formalism 
suffers from the neglect of the effects of the ionized medium. 
Therefore , in order to improve our understanding of the nature of 
these solar bursts it would be of considerable importance to re- 
lax the ultrarelativistic approximation used in previous treatments 
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and to evaluate synchrotron emission spectra from electrons of 
arbitrary energies in the presence of an ionized medium. 

The general problem of the radiation from charged particles 
moving in a magnetoplasma was treated by Eidman (1958) , Liemohn 
(1965) and Mansfield (1967). These treatments provide explicit 
expressions for the spectral and angular distribution of the 
power radiated into both ordinary and extraordinary modes by a 
single electron moving in a helical orbit in a cold, collision- 
less plasma permeated by a static uniform magnetic field. The 
resulting spectra were numerically analyzed by Liemohn (1965) at 
frequencies in the vicinity of the plasma and gyro-frequencies of 
the ambient electrons for which the indices of refraction of either 
the ordinary or extraordinary modes are real and greater than unity. 
Although part of the solar radio emission, especially in the micro- 
wave band, may be produced at these frequencies, in the present 
paper we shall only investigate the radiation at frequencies above 
the plasma frequency where the indices of refraction of both 
ordinary and extraordinary modes are real and less than unity. 

For such frequencies the phase velocity of light is greater than 
c, and, as pointed out above, this results in the suppression of 
the low frequencies of synchrotron emission. In order to assess 
the importance of this effect for intermediate energy electrons, 
an expression is derived for the frequency spectrum, integrated 
over all angles, of the emission from an electron of arbitrary 
energy moving in a circular orbit in a homogeneous and isotropic 
electron plasma. This derivation is based on the more general 
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treatments mentioned above, as well as on an independent de- 
rivation appropriate for circular motion and frequencies for 
which both indices of refraction are isotropic, real, and less 
than unity. The resultant emission spectra are evaluated nu- 
merically for various electron energies and plasma parameters 
and it is shown that for large plasma frequencies there is a 
significant suppression of emission for electrons of all energies. 
Finally, by comparing these spectra with those obtained from the 
ultrarelativistic approximation mentioned above, the validity of 
the high-energy formulas used in the previous treatments of the 
effects of the medium is investigated at various electron energies 
which could be of importance for solar radio emission. 



Radiation From Arbitrary Energy Electrons 


The angular and frequency distribution of the electro- 
magnetic emission from an electron moving with an arbitrary 

velocity v in a circular orbit in a homogeneous ionized medium 

£ 

permeated by a static uniform magnetic field u is given by (Liemohn , 


1965) : 
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The subscripts (+) or (-) indicate emission into the ordinary 
or extraordinary modes respectively; g is the ratio of v to c; 

Y is the electron Lorentz factor; 6 is the angle between ^ and 
the radius vector from the electron’s guiding center to the point 
of observation ; is the cyclotron frequency of the ambient elec- 
trons ; D +(v , 6) is the index of refraction of the ionized medium 
and is in general frequency dependent and anisotropic with respect 
to the direction of the magnetic field; and Jg is a Bessel function 
of order s. The polarization coefficients oH and 

(Liemohn , 1965) are defined in terms of the components of the 


electric vector of the radiation field 
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where E and Eg are the transverse components of E and Ejj. is the 
component of E along the 9 direction defined above. 

The index of refraction and the polarization coefficients 
can be determined from the properties of the ambient ionized 
medium. For a cold collisionless electron plasma these quantities 
are (Ginzburg , 1961; Liemohn , 1965). 


i 2_P^fP^Zfl}- 

" ~ / * ±rp‘ t sm‘ , & * ‘if x (P l -F l ) X co^^' l3 --lFHf-F l ) - F l Sin& 
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The dimensionless quantities F and P are defined in terms of the 
cyclotron and plasma frequencies 

F * W*>e j P « >> p / v e 


where Vb and are given in terms of the magnetic field B and 
the electron density fl^by 

Vq 2 1 1 2J? efr/MC ; \) p = 

The indices of refraction , #7+ and H - , have several cutoffs (*»=o) 
and resonances . For FyP»*Tl/^^ > however, ft + is real and 
less than 1, and for F > ^ 
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is also real and less than 1. Assuming that and limiting 

our discussion to frequencies greater than ( 


we see that for both indices of refraction are real and there- 

fore the total emission is the sum of the emissions in the ordi- 
nary and extraordinary modes given by equation (1). Since 
we can neglect second order terms in 1/F and 1/P, and therefore 
the index of refraction and the polarization coefficients given 


by equations (2) , (3) and (4) reduce to 

** - / - 


(£) 


= O 


Using equations (5) , (6) and (7) , the total emission is obtained 
by adding the contributions from the two modes in equation (1). 


This results in 


Jiuv 


~ £2 T 


2 

Jco^O £ (</>»**&)] f(*~ J m ) 


In order to substantiate the validity of the limiting process used 
to derive equation (8) , we provide an independent derivation of 
this equat ion j appropriate for the physical conditions mentioned 
above. This is done in Appendix I, where we show that equation 
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(8) can indeed be derived directly from Maxwell's equations for 

an electron of an arbitrary velocity moving in a circular orbit 

in a static uniform magnetic field immersed in an isot ropic 

electron plasma with index of refraction given by equation (5). 

Since the radiation has azimuthal symmetry with respect to 
— * 

the direction of B, the total radiation, I CP) , integrated over all 
angles is given by 

X(v) - 2Trf ^ S'* & Of) 

This integral can be evaluated by the same methods as those used 
by Schott (1912) for the evaluation of a similar integral for 
radiation in vacuum. The presence of the ionized medium, however, 
introduces some complications and therefore in Appendix II we pro- 
vide an outline of the integration process suitable for the pre- 
sent problem. The resulting integrated emission is given by 

f\(***) J *] < !(»-*?») to) 

For fl=l equation (10) directly reduces to the emission spectrum 
in vacuum given by Schott (1912) and Landau and Lif schitz (1962). 

In order to compare equation (10) with previous treatments of the 
effects of the medium on synchrotron emission we introduce the 
emissivity function, F(»>/*) c ) , given by 
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where 


iL - hJj - ii. 

Vc. ~ 3 / 3 


in terms of which equation (10) can be written as 




(/ 2 ) 

fa) 


For certain limiting conditions equation (11) can be considerably 
simplified by using asymptotic forms of the Bessel functions. 


In Appendix III we show that: 

a) If s '*;»• equation (11) reduces to 


f(±) -[* Q jl r>-i */>,) 
lv < ! U *» >W 777 f J 


b) If S»l and V, » I (but S/ tf not necessarily larger 


than 1) equation (11) also reduces to 
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From equation (17) we see that, depending on the value of y 

condition (a) can be satisfied for small as well as for large 
values of s and y. For example, for =*.S/> 3 for all 

values of s and y of interest. On the other hand, condition (b) 
will be satisfied in general only by ultrarelativistic electrons. 
Indeed, equation (15) is identical to the emissivity function 
used in the previous treatments of the influence of the ionized 
medium ( Ginzburg and Syrovatskii , 1964, 1965; Ramaty and Lingenf elter , 
1967.) The validity of this approximation, however, can be best 
assessed by the numerical evaluations which are discussed in the 


next section. 
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Numerical Results 

The emissivity function F (»/*) can be evaluated by using 
equation (11), or if either conditions (a) or (b) are satisfied, 
F(u/v c ) can be directly obtained from equations (15) or (16) 
respectively. The use of these asymptotic equations is necessary 
in certain cases owing to the difficulties involved in the nu- 
merical evaluation of Bessel functions of large orders and argu- 
ments. By simultaneously evaluating equations (11) and (15) we 
find that the values of F ( u /‘0 given by these equations do not 
differ by more than 20% if 3 • Similarly, we also find 

that equation (15) is accurate to within 20% if ^ > 2 - 5 - Using 
these transitional conditions we have evaluated as a 

function of v/>4 for various values of ^ and V»a 

In figure 1 F (»/v c ) is plotted for l)p = o (emission in vacuum.) 
The circles correspond to the first 10 harmonics (s=l,...,10) and 
it can be seen that as the electron energy increases the harmonics 
cluster together and the emission spectrum. can be approximated by 
a smooth distribution. It can also be seen that the emissivity at 


high frequencies increases with increasing energy and that for 
Lorentz factors greater than about 3 , F ( Vv c ) can be well approxi- 
mated by the ultrarelativistic spectrum C^"* 00 ) calculated from 
equation (15) with H = l. 


In figure (2) F (v/» 4 > is plotted as a function of for 


V v 6 = 6 ' 

of equation 


(The assumption that V>S made in the derivation 
(8) is quite well satisfied in this case.) As in the 


previous ultrarelativistic treatments of the effect of the medium 
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(e.g., Ramaty and Lingenf elter , 1967, Figure 1), there is a sig- 
nificant suppression of the emission at low frequencies. Moreover, 
for low values of y, the total radiated power is also strongly 
suppressed. This total suppression is the combined result of the 
low frequency cutoff due to the influence of the ionized medium 
and the small emissivity of the mildly relativistic electrons at 
high frequencies. 

In order to assess the deviations of the ultrarelativistic 
approximation mentioned above from the present results we have 
plotted in figure (3) the ratio, Fjjr/F , of the spectra obtained 
from equation (15) to those shown in figure (2) . We can see that 
the ultrarelativistic approximation is quite adequate at large values 
of y since for all values of where there is significant emission, 

F ur /F does not deviate significantly from unity. At lower values of 
y , however, the ultrarelativistic approximation appreciably over- 
estimates the intensity of the emission at all values of »/i> c • 

Because of this, for sufficiently large plasma frequencies, the 
cutoffs in the emission spectra of intermediate energy electrons 
due to the influence of the medium will be even sharper than the 
low frequency suppressions obtained in the previous ultrarelativistic 
treatments of this effect. As pointed out in the introduction, this 
may have interesting implications on Type IV solar radio bursts. 

In figures (4) and (5) we plotted for V* of 3 and 

1.5 respectively. It can be seen that the influence of the medium 
diminishes with decreasing plasma frequency. The approximations 
based on the assumption that V> v e , mentioned above, probably 
break down for the spectra shown in figure 5. However, as can be 
by comparing figures 1 and 5, for =1.5 the influence of 


seen 
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the medium becomes quite small and at frequencies where both modes 
propagate freely the emission is much like that in vacuum. 
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Appendix I 

We consider the radiation produced by an electron moving in 

— » 

a circular orbit in a static and uniform magnetic field B immersed 
in a homogeneous and isotropic medium characterized by a dielectric 
constant £ and a magnetic permeability equal to that of free space. 
The electric and magnetic vectors of the radiation field satisfy 
Maxwell's equations: 


-• 

V-E 

- HTf / e 

(11) 

7-H 

* 0 

(12) 

Vxi$ 

c J c 

0-5) 


+ s 0 

+ c 

(«) 


where 0 and J are the charge and current densities associated with 
the electron's motion. A wave equation for the magnetic vector H 
can be obtained by taking the curl of equation (13) and by using 
equations (12) and (14): 

f ”7 (15) 


Since the motion of the electron is periodic, with gyroperiod T, 
both H and 7 can be decomposed into Fourier series. 




H (?,*)« Z e ' T 

s 
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(Ji) 
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the coefficients of which satisfy Poisson’s equation: 

+ V* H s = - f Cl8) 


where 


U . nrS'fl 
** " ”Tt“ 


(i ‘ i) 


By using the inverse transform of (17) , the solution of (18) can 
be written as 


, . t cV s l?-M 

* f- [it e ^ [iv ?»!(?') - e i ft ^ lo) 


Since the current density can be written as 

* e vft) £ ( ? - ? «) 


(m) 


where (t) and V (tj are the electron’s instantaneous position 
and velocity, equations (110) reduces to 

O 

-*■ 

where the gradient is taken with respect to ^ . 


We now define a coordinate system such that its origin is at the 
electron’s guiding center and that the z axis is parallel to the 
direction of the static field B. Since the radiation is symmetric 
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with respect to B, we can define the azimuthal orientation of the 
system so that the yz plane contains the point of observation 
Let Ro be the distance from the origin to the point of 

A 

observation and I? a unit vector along this direction. If R 0 is 
much larger than both the wavelength of the radiation and the 
gyroradius fg , | r ~ I can be approximated by 

A 

I r - r e | « R 0 -*?e • I? 


Using this relation and neglecting • I* in the denominator of 
the term in square brackets, equation (112) can be written as 

(I >3) 




where 


and where 


t -U 

s 


\Cr) 


R (?) - & — 


is given by 

[ S*o 


T 


[it 5 ft) 


c [ 


O) 


In the coordinate system defined above, the z component of A g 
vanishes and its x and y components are given by 

v[s'f-(ep si »»& 4 inV f] 






(J»s) 
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ft 


- ei t 
s t ire R. 


(xi 6 ) 

o 

where Tj is the magnitude of V* , is the ratio of V to c and 9 
is the angle between W and B. Using some well known integrals 
from the theory of Bessel functions, (115) and (116) reduce to 


ftsx « Hr t 3.' ( sp* •*»«!) 


cR 
2e 


9 sh = 


j s (sp|5«.x&) 


(xr?) 
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In an isotropic dielectric medium the electric and magnetic 
vectors of plane electromagnetic waves and the propagation vector, 

A 

^ , are mutually orthogonal and related by 


e -£*->• 


Cm) 


As discussed in the main text, however, this relation is valid 
only for sufficiently high frequencies where the anisotropy in- 
troduced by the static field becomes negligible. Using (119), 


the pointing vector £? can be written as 

s = 


clHl* ^ 


*nr\fe. 


(p-) 


and therefore, the instananeous power radiated into the solid 
angle element is is given by 


Mir vt 


(1*0 
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The average power radiated into the s harmonic is obtained by 
using (16) and by averaging (121) over one gyroperiod T. This 
results in 


JU = 








(l22.) 


Making use of equations (113) , (117) and (118) , equation (122) 
can finally be written as 


ai 4 _ 

' 




£E£Z S 

cT‘ 


» 

\fe 


cot*e fjspfi s-«») + ft £ (sf.'li «» »)] 


Since is equal to the index of refraction n, and the gyro- 

period T is given by 



equation (123) is to equation (8) given in the main text 
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Appendix II 


We now provide an outline of the evaluation of the in- 
tegral in equation (9). Using the identities ( Schott , 1912; 
Watson , 1966) 

IT 

Jj <■)>> = 4: J J o (ixsiof) 3.s<f J<f (ffl) 


(*) = £ j J. (2*Si« f) ( coslf--^)on2sf Jf (Jz) 

0 

the term in square brackets in equation (8) reduces to 

'T 

~ \ J 0 CnZS'f [ pV-l . 


Using this expression and the integral ( Gradshteyn and Ryzhik, 


1965) 


V t - *** 'j. 


equation (9) , with ji nsi given by equation (8) , reduces to 

w») = - seia 1 . 
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Making use of the identities ( Schott , 1912; Schwinger , 1949) 

jJ. Li) - - k- \ Sl '« C»2sf ttnf 


L(K)J» - if *"!?*!?.*> 
* i 



equation (II 13) becomes equivalent to equation (10) . 
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Appendix III 


We shall first derive equation (14) by using an asymptotic 
formula from the theory of Bessel functions ( Watson , 1966) , 


3/4 


namely that for S C * ~ 


J 2 , ( 2 san) « -1- f J 

2S ' {Ws (i-ti ri ) L -* 


2.S 


{IfFs 


By differentiating and integrating this expression and by 
using the identity 



, which can be derived by direct differentiation, and the con- 
t 

dition J I , mentioned above, we find that 

J ' ttfi H [-EL. f * (if i) 


and 
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Using equations (III1) and (III2) , equation (11) reduces to 
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Finally , making use of equation (17), we obtain that 


F (*«) *Ji 




As can be seen, this expression is identical to equation (14) 
in the main text. 

We now consider the derivation of equation (15). Making use 
of the asymptotic formula (Watson, 1966) 




a(s-y) 

ax 



3*" 1 


J 



which is valid for large s and x, and x< JJ, we find that 




^ [ % ( l 
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(fffij) 


By integrating and differentiating equation (III6) we directly 


obtain 


J I** (*SX) « 






V hL 


and 






where £ is given by equation (16). 
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In the process of differentiating equation (III6) , we made use 
of the formula ( Schwinger, 1949) 

It K.^60 + vt^(x) = -V *■**(*) 


Substituting (III7) and (III8) into equation (11) and again 


making use 

of the condition 

, we obtain 




to 




- j 

J C® *0 





Using the identity ( Schwinger , 1949) 


equation (I I 19) can be reduced to 

FUJ * -*(^) j i , >r-»* (H'°) 

Wv* 

which is identical to equation (15) given in the main text. 

As pointed out above, equation (III4) is valid for all values 
of s and ^as long asS/j\>!is satisfied. If, however, ^r f 
also is much larger than unity, equation (III4) should reduce 
to an appropriate asympotic form of equation (III10) . Indeed, 
if ^ » | , the expression in square brackets in equation 

(III4) can be expanded in powers of ^ : 

2 

- To> 


, 2 . 

* e. » I - * 
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Using this expression, and assuming that H I , equation 
(III4) reduces to 



By making use of an asymptotic expression given by Ginzburg 
and Syrovatskii (1964), namely that for 

x [ (l) 

* 


we see that equations (III4) and (III10) have indeed the 
same asymptotic forms for Jf ( »| and S>^. f 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


FIGURE CAPTIONS 


The emissivity function '(«/*) in vacuum for electrons 
of various Lorentz factors. The circles denote the 
frequencies corresponding to the first 10 harmonics. 

The actual emission consists of discrete spikes at the 
individual harmonics and can be obtained from F (•/Hr) 
by using equation (13) . 

The emissivity function F (*/v c ) in an ionized medium 
for \)p and for electrons of various Lorentz factors. 

The actual emission consists of discrete spikes at the 
individual harmonics and can be obtained from F (o/vj 
by using equation (13) . 


The ratio between the values of F (v/i> c ) calculated by 
using the ultrarelativistic approximation to those shown 
in figure 2. F^ corresponds to the values of F (v/V c ) 
obtained from equation (15) for l ^5 = 6 . 


The emissivity function F 6>/v c ) in an ionized medium for 
V Iv and for electrons of various Lorentz factors. The 

r ft 

actual emission consists of discrete spikes at the in- 
dividual harmonics and can be obtained from F (»M) 
by using equation (13) . 


: The emissivity function F^u/v c ) in an ionized medium for 

_ j.5 and for electrons of various Lorentz factors. 
The actual emission consists of discrete spikes at the 
individual harmonics and can be obtained from F (v/v e ) 
by using equation (13) . 
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